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ABSTRACT 

The anomalous lack of large angle temperature correlations has been a surprising feature of 
the CMB since first observed by COBE-DMR and subsequently confirmed and strengthened 
by WMAP. This anomaly may point to the need for modifications of the standard model of 
cosmology or may show that our Universe is a rare statistical fluctuation within that model. 
Further observations of the temperature auto-correlation function will not elucidate the issue; 
sufficiently high precision statistical observations already exist. Instead, alternative probes 
are required. In this work we explore the expectations for forthcoming polarisation observa- 
tions. We define a prescription to test the hypothesis that the large-angle CMB temperature 
perturbations in our Universe represent a rare statistical fluctuation within the standard cos- 
mological model. These tests are based on the temperature- Q Stokes parameter correlation. 
Unfortunately these tests cannot be expected to be definitive. However, we do show that if 
this TQ-correlation is observed to be sufficiently large over an appropriately chosen angular 
range, then the hypothesis can be rejected at a high confidence level. We quantify these state- 
ments and optimise the statistics we have constructed to apply to the anticipated data. We find 
that we can construct a statistic that has a 25 per cent chance of excluding the hypothesis that 
we live in a rare realisation of ACDM at the 99.9 per cent confidence level. 

Key words: cosmology: cosmic background radiation - cosmology: large-scale structure of 
Universe 



1 INTRODUCTION 

In the two decades since the Cosmic Background Explorer (COBE) 
first detected the primordial fluctuations in the cosmic microwave 
background (CMB) temperature jWright et al.|1992) , and perhaps 
even more so in the last decade over which the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) has provided ever more ac- 
curate full-sky maps of those fluctuations (see |Komatsu et al.|201 1| 
for example), the CMB has become a keystone in the remarkable 
transition of cosmology from a qualitative to a precision science. 

An important element of the role of the CMB in precision cos- 
mology has been that the canonical theory of cosmology, inflation- 
ary Lambda Cold Dark Matter (ACDM), makes clear predictions 
for the statistical properties of the spherical harmonic coefficients 
of the temperature fluctuations, 

ae^ = J Y;^{e,<l})T{9,<j,)d{cose)d<j>, (1) 
which are predicted to be statistically isotropic realisations of inde- 
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pendent Gaussian random variables of zero mean and with variance 
Ce depending only on £, 

I 

= ^ E (2) 

In modem discussions of the CMB the two-point angular power 
spectrum, embodied in these Ce, plays a central role and is the 
source of the remarkable precision of the cosmological parameters 
( |Komatsu etal.|2011t . 

Before the COBE era, it was the two-point angular coiTelation 
function of the fluctuations, 

C{9) = r(ni)r(n2)ifti.ft2=cos9, (3) 

rather than the angular power spectrum that was of primary interest 
to astronomers. Statistically the two-point angular correlation func- 
tion is an ensemble average but, in practice, this must be replaced 
by an average over pairs of points separated by an angle 9, as de- 
noted by the bar over the expression. In fact the COBE Differen- 
tial Microwave Radiometer (COBE-DMR) did report C{9), though 
only in their final, four-year paper ( Bermett et al.|1993) . When ex- 
tracted from a full sky map both C{9) and the Ce contain the same 
information, albeit in different forms. The same is true between a 
function and its Fourier transform; signals are typically most easily 
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seen in one or the other forms but not both. In the case of the CMB, 
C{9) and the Ct are related by a Legendre transform. The Ce most 
easily show the small angular scale, microphysics at last scattering; 
whereas C{0) most easily shows the large angular scale behaviour. 

As observed by the COBE-DMR, C{6) had an unexpected 
property - it was consistent with zero for angular separations be- 
tween approximately 60° and 160°. This was duly noted at the time 
but mainly remembered today as a low quadrupole. The WMAP 
team confirmed the COBE-DMR observation of a lack of large- 
angle correlation with significantly smaller error bars. In their ini- 
tial, one-year release, [Spergel et al.| ( j2003j phrased the anomaly in 
terms of a statistic 



•5*1/2 = 



1/2 



(4) 



In the best fit ACDM model the expected value of this statistic is 
approximately 50,000 {\iK)^, whereas the observed value is ap- 
proximately 8500 (viK)^ on the full sky, for example from the 
WMAP Independent Linear Combination (ILC) map[^ with a p- 
value of approximately 0.05. Even more striking is that if one con- 
siders only the part of the sky outside a conservative galaxy cut, 
then 51/2 ^ 1000 - 1150 (^K)^ and is only ~ 1300 (^K)^ in 
each of the V and W frequency bands, which are expected to be 
dominated by the CMB signal. The cut-sky S1/2 has a p-value less 
than about 2.5 x 10^*, depending on the precise map ( Copi et al. 
[2009l l. 

We have argued that such absence of two-point angular corre- 
lation is unlikely to result solely from a small quadrupole, or even a 
small quadrupole and octopole, and that it instead requires a "con- 
spiracy" among the first several multipoles ^Copi et al.|2009^ . Such 
covariance among the C'e is likely contrary to the fundamental pre- 
diction of the canonical cosmological model that the underly- 
ing the Ci are Gaussian random independent variables with vari- 
ances depending only on £. This is one of a number of large scale 
anomalies that suggest modifications of the standard model are re- 
quired on large angular scales (seel Bennett et al.|2011[[Copi et al.| 
|2010 . and references therein for further details). 

One possible explanation of the absence of large angle corre- 
lations among the CMB temperature fluctuations is that it is merely 
a statistical fluke. In this paper we explore the consequences of this 
hypothesis. In particular, since the fluctuations in the CMB tem- 
perature and in its polarisation arise largely from the same source 
- the gravitational potential - one might have hoped that a small 
temperature-temperature (TT) correlation function on large an- 
gular scales would predict a similarly small cross-correlation be- 
tween CMB temperature and CMB Stokes parameter Q, or in the 
polarisation-polarisation (QQ) correlation. 

Unfortunately, as we shall see, the connection between tem- 
perature and polarisation fluctuations is too weak for a general 
definitive test of the origin of the vanishing correlation function. 
However, we do find that if the 51/2 is small because of a statistical 
fluke within ACDM cosmology, then the cross correlation between 
temperature and polarisation is unlikely to be large on large angu- 
lar scales. Therefore were we to infer a large value of this cross- 
correlation from future data, that would be evidence against a sta- 
tistical fluke as an explanation of the vanishing TT correlation. 

In this paper we provide a prescription to follow in order to 
test this hypothesis. In Section [2] we describe the construction of 
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an ensemble of realisations of ACDM that are constrained to re- 
semble our observed universe in the properties of their TT an- 
gular power spectrum and full-sky and cut-sky two-point angular 
correlation functions. In Section[3]we construct statistics that, like 
5*1/2 for the TT correlations, can be used to quantify the small- 
ness of the correlations between temperature and polarisation fluc- 
tuations. Section |4] contains a discussion of the results of applying 
the temperature-polarisation cross-correlation statistics to the en- 
semble of constrained realisations and looks forward to what might 
be learned by applying them to future polarisation data. Given the 
constructed realisations the value of the statistics are fixed and the 
optimal application of them to polarisation data can be determined 
as discussed in this section. This optimisation is independent of po- 
larisation observations. Finally, Section|5]contains the conclusions. 



2 CONSTRAINED REALISATIONS 

To study the signature of the lack of large angle correlations in the 
WMAP temperature data on upcoming polarisation measurements, 
such as from Planck, we require realisations of ACDM consistent 
with these large angle results. For this purpose we have generated 
300, 000 such realisations as follows. 

(i) The Cf ^ at low-^ are well measured. Assuming Gaussian 
statistical errors we generate random Cj'^ from Gaussian distribu- 
tions centred on the WMAP-reported values. This produces a power 
spectrum consistent with that reported by WMAP. It is important to 
stress that this is a power spectrum consistent with the observa- 
tion of our particular realisation of the Universe as measured by 
WMAP, not a general realisation of the best fit ACDM model. For 
this reason cosmic variance is not relevant nor do we generate Ct 
from a distribution. It is also true that on a partial sky the Ci are 
slightly correlated. To correct for this we actually use the Fisher 
matrix from the WMAP likelihood code without the contribution 
from cosmic variance in drawing these Cj"^ . In practice this is a 
small correction but has been included for completeness. 

(ii) Given the Cj'^ from the previous part and assuming Gaus- 
sian random phases, a single sky realisation, the aj^, is generated. 
The resulting sky realisation has a 5i/2 consistent with the small 
value in the full sky WMAP ILC map. 

(iii) To further be consistent with WMAP observations the 5i/2 
on the cut sky must also be small. For our realisations we require 
5i/2 1292.6 the value from the WMAP seven-year, 
KQ75y7 masked ILC map. This 5i/2 value is calculated for a re- 
alisation by first constructing a map at NsiDE=64 from the aj^ 
generated above. The pseudo-C^ are extracted from this map based 
on the region outside the KQ75y7 mask using SpICE ( jChon et ^ 
2004 1. Finally 5i/2 is calculated using these Ci up to l-max = 100. 



(iv) For temperature realisations that satisfy the cut sky con- 
straint we also generate realisations of the af^. See Appendix [a| 
for a review of the process. 

The construction of a set of constrained realisations is the basis for 
the prescription we are describing. It will next provide predictions 
for the expectations from the observations of the CMB polarisation. 



3 STATISTICS 

For the temperature auto-correlation the 5i/2 statistic was defined 
a posteori ( |Spergel et al.|2003| l to be 
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Figure 1. Example histogram of the 5-^*3 statistic jSj for constrained (solid, 
black line) and ACDM (dashed, red line) realisations. We note that the 
constrained realisations are more sharply peaked at low 5^*^ than ACDM 
which, though peaked at approximately the same value, has a long tail. The 
dashed, vertical lines represent the values of S'"^'^ for which 99 per cent 
and 99.9 per cent, respectively, of the constrained realisations have smaller 
values. 



Figure 2. Example histogram of the 5^*3 statistic jo} for constrained (solid, 
black line) and ACDM (dashed, red line) realisations. We note that the con- 
strained realisations are more sharply peaked near zero than ACDM which, 
though also peaked near zero, has a long tail particularly to large, positive 
values. The dashed, vertical lines represent the values of s'^'^ for which 99 
per cent and 99.9 per cent, respectively, of the constrained realisations have 
smaller values. 



S,/2= J [C^^(e)]'d(cose). (5) 

Inspired by this we define a comparable statistic for (7"^'^(S), the 
two-point angular correlation function between fluctuations in the 
temperature and the Stokes parameter Q. 

Observable properties of photons can be characterised by the 
Stokes parameters. For the CMB the relevant quantities are the in- 
tensity, conventionally represented by the temperature, T, and the 
linear polarisation given by the Q and U parameters. For the CMB 
the circular polarisation, represented by the V Stokes parameter, is 
expected to be zero and not considered further. When working in 
real space the natural correlations to construct are amongst these 
observables, T, Q, and U. These correlations are constructed such 
that they only depend on the angular separation between pairs of 
point on the sky and are thus rotationally invariant despite the fact 
that the definition of Q and U depend on the choice of coordi- 
nate axes. When working in harmonic space, it is natural to de- 
compose the polarisation into "gradient" and "curl" modes alter- 
natively called E and B modes, which are similarly rotationally 
invariant quantities. These latter names will be used throughout. 
Thus in real space we will work with the TQ two-point angular 
correlation function, C'^''^{9), which may be written in terms of 
the two-point angular power spectrum coefficients, Cj^ . Details 
of these representations can be found in standard references i Zal-| 
[darriaga & Seljak|1997{[Kamionkowski et al.|1997[ for example). 

3.1 S'^'^ Statistic 

In the case of polarisation a priori the optimal range over which to 
integrate the correlation function is unknown and will be explored 
below so we define the general statistic 

/•cos 6i 

s^'^ie^,e2) = / [c^«(e)]'d(cose). (6) 

•J cos 02 



As with Si/2 we may calculate this easily in terms of the power 
spectrum coefficients, Cf^. Using 

oo I— — 

we may show that 

s^''{e,,e,) = Y.cril^{e,,9,)cr, (8) 

t,e' 

where ij^ {61,92) are components of a known matrix calculated 
in Appendix |b] A histogram of the S^*^ statistic for a particular 
choice of 9i and 6*2 is shown in Fig. [T| comparing the constrained 
realisations to ACDM. 



3.2 3^*5 Statistic 

Motivated solely by its simplicity and ease of computation we also 
define a new statistic which is linear, rather than quadratic, in the 
TQ correlation function 

/"COS 61 

s^'^{9^,92)= / C'^'^{9)A{cos9). (9) 

J cos 0^ 

As with S'^'^{9i, 02 ) we may calculate this easily in terms of the 



s'^'^{9i,92) ^^Cf if '^(91,92), (10) 

e=2 

where if^ {9i , 6^2) are the components of a known vector calcu- 
lated in Appendix|B| A histogram of the s^'' statistic for a particu- 
lar choice of 9i and 92 is shown in Fig.|2]comparing the constrained 
realisations to ACDM. 
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4 RESULTS 

The S^''^ and s"^*^ statistics defined above iiave been calculated 
for the constrained realisations discussed in Section |2] and for a 
comparable number of realisations of ACDM. In both cases these 
have been calculated from maps produced at NsiDE=64. Data from 
the temperature map outside the KQ75y7 mask and from the seven- 
year polarisation analysis mask, both provided by WMAP, has been 
used. As shown in Figs. [T| and |2] the constrained and ACDM real- 
isations have most likely values for the S"^*^ and s^*^ statistics at 
nearly the same value. However, as we also see ACDM predicts 
much broader distributions for the two statistics. In particular there 
is a significant probability of producing values larger than the con- 
strained realisations. This provides a means of testing the hypothe- 
sis that our Universe is just a rare realisation of ACDM. This results 
in a simple but not definitive test. 

Consider the case of the S'^^ statistic as represented in 
Fig. [T| For the constrained realisations 99 per cent of them 
have S'^'3(48°,120°) ^ 1.403 (^K)* and 99.9 per cent have 
S'^'^ (48°, 120°) < 2.195 (^K)■*. Unconstrained ACDM (with 
the best fit values of cosmological parameters) randomly gener- 
ates realisations with values larger than these 38.6 per cent and 
25.6 per cent of the time, respectively. If observations of the polar- 
isation show our Universe to have a 5^*2(48°, 120°) value larger 
than these values then we can reject the random ACDM realisation 
hypothesis at the appropriate confidence level. Alternatively, if the 
5'^<3(48°^ 120°) value is smaller then no definitive statement can 
be made; the polarisation fluctuations would be consistent with the 
hypothesis that we live in a rare ACDM realisation but do nothing 
to advance that hypothesis. This is the main point of the paper. 

Similar statements may be made about the s^'' statistic shown 
in Fig. [2] It provides similar information as S"^^ . 

It remains to determine the optimal ranges over which to eval- 
uate the statistics. We define optimal to mean the maximum dis- 
criminatory power between the distribution of the statistic in the 
constrained versus ACDM realisations. For a given per cent cutoff 
from the constrained realisations we wish to find the angle range 
[9i , 62] that has the maximum fraction of ACDM realisations above 
this value. 

The results of such a study are shown in Figs.j3]and|4]for the 
S^'^ {61 , 62) statistic and in Figs.[5]and|6]for the s'^'-^6i , O2) statis- 
tic. The statistics are non-zero only up to the diagonal 9i = 62, 
shown as the dashed, grey line in the figures, but not along it. For 
this reason the contours are truncated at the diagonal. They have 
been made symmetric in 9i and 62 (by taking and Is"^*^!) so 

the results are shown as identical when reflected through the diag- 
onal. The optimal ranges and fractions of ACDM realisations are 
listed in Table [T| We note that the optimal surfaces represented by 
the contours seen in the figures are relatively broad, at least in one 
direction. Due to this the values of 9i and 62 in a neighbourhood of 
those listed in Table[T|can be employed with equal efficacy. 



5 CONCLUSIONS 

The absence of two-point angular correlations on large angular 
scales in the CMB temperature data is, by now, well established. 
It was first measured by COBE-DMR, but became more significant 
in the WMAP temperature maps. This absence is difficult to accom- 
modate within the standard cosmological model, especially since it 
seems to imply co- variance among low-^ multipoles of the CMB. 
A simple explanation that has been proffered is that we just live in a 




01 (deg) 



Figure 3. Contours for the fraction of ACDM realisations above the 99 per 
cent value of the constrained realisations from the statistic jS). In the 
optimal case 38.8 per cent of the ACDM realisations have a larger value. In 
this figure the results around not defined along the diagonal (grey, dashed 
line) and have been made symmetric about it by taking the absolute value 
of the statistic. 
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Figure 4. The same as Fig.[3]now for the 99.9 per cent case. Here in the 
optimal case 25.6 per cent of the ACDM realisations have a larger value 
and the full histograms are shown in Fig.[T] 
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Figure 5. Similar to Fig. [s] but now for tlie s-^*? statistic j9}. Here in the 
optimal case 38.5 per cent of the ACDM realisations have a larger value at 
the 99 per cent level. 



Figure 6. The same as Fig.|5]now for the 99.9 per cent case. Here in the 
optimal case 26.4 per cent of the ACDM realisations have a larger value 
and the full histograms are shown in Fig.|2] 



Table 1. Optimal angle ranges for the S^<5 statistic (61 and s'^'^ statis- 
tic |9j. The optimal ranges are determined by finding the maximum fraction 
of ACDM realisations with the appropriate statistic above the 99 per cent 
or 99.9 per cent level of the constrained realisations. The histograms for the 
optimal 99.9 per cent ranges are shown in Figs.[T]and|2] The full contours 
are shown in Figs.[3]j6] 



Statistic 


C.L. 


01 


02 


Fraction 




(per cent) 


(degree) 


(degree) 


(per cent) 


STQ 


99 


46 


122 


38.8 




99.9 


48 


120 


25.6 




99 


20 


177 


38.5 




99.9 


23 


174 


26.4 



rare realisation of ACDM that happens to have a lack of large angle 
TT correlations. If so, one might hope that constraining ACDM re- 
alisations to have low TT correlations at large angles, would have 
observable consequences for other correlation functions, such as 
TQ. These would be the basis for an observational test of this sta- 
tistical fluke hypothesis. 

In this paper we have discussed a prescription to follow in 
order to test this hypothesis of our Universe being a statistical fluke. 
The prescription may be simply stated as follows. 

(i) Construct realisations of our Universe consistent with the ob- 
served temperature fluctuations. This means construct sets of aj^ 
and af^ consistent with the observed Cj"^ and full and cut sky 
5*1/2 as discussed in Section[2| 

(ii) Apply the 5'"^'^ and statistics as defined in Section|3]to 
these constrained realisations. 

(iii) Also apply the S*"^*^ and s'^'^ statistics to a comparable 
number of best fit ACDM realisations and use these to find the 
optimal ranges [6i, 82] for each statistic. Optimal here means that 



the maximum fraction of ACDM realisations fall above the value 
at some confidence level in the constrained realisation, for example 
the 99 per cent or 99.9 per cent level. 

(iv) Given the optimal ranges from the previous step now ap- 
ply these particular cases to the observed polarisation signal. If the 
observations produce values for these statistics larger than that ex- 
pected from the constrained realisations, then the statistical fluke 
hypothesis can be rejected at the appropriate confidence level. Al- 
ternatively if the values are larger, then the hypothesis remains con- 
sistent but unproven. 

We further note that the optimisation in this prescription is inde- 
pendent of the polarisation observations, or, in fact, whether the 
polarisation has been observed or not. 

Our work is, in spirit, related to |Dvorkin et al.|p008 t Whereas 
they consider observables in the polarisation signal for "models" of 
three dimensional primordial power modulation that might explain 
the breaking of statistical isotropy in the temperature field, we pre- 
dict the polarisation statistics starting directly from realisations of 
ACDM models that are constrained to show the suppressed TT 
correlation at large angular scales. 

In the work reported here we have generated realisations and 
performed the optimisation based on the WMAP seven-year data 
release. The prescription could be applied to the WMAP nine- 
year data release and the results are not expected to differ signif- 
icantly. We have also said nothing about applying the statistics to 
the WMAP reported polarisation observations. Unfortunately the 
signal-to-noise in the polarisation observations is not yet sufficient 
to make meaningful statements. To see thi s, using the WMAP nine- 
year reported Cf^ ( Hinshaw et al. |2012| we find for the optimal 
ranges 



qTQ 

^WMAP 
and 



(48°, 120° 



(1.0 ±0.8) (^K)" 



(11) 
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TQ 

^WMAP 



(23°, 174°) = (0.8 ±0.8) (^K)' 



(12) 



Here the error bars are crude estimates assuming the reported 
are statistically independent and the noise is Gaussian. These as- 
sumptions are not justified and a more careful assessment could be 
performed using the Fisher matrix. However, given the large esti- 
mated errors such an assessment is not warranted. 

We have shown that the prescription described in this work is 
not a definitive test of the statistical fluke hypothesis for our Uni- 
verse. Nevertheless, by carefully optimising the statistical measure 
of large angle TQ correlations, we were able to demonstrate that 
once good TQ correlation data is available there is a reasonable 
probability (over 25 per cent) to reject the statistical fluke hypoth- 
esis at the 99.9 per cent confidence level. WMAP data is not up to 
this task, however Planck data should be. 
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APPENDIX A: POLARIZATION GAUSSIAN RANDOM 
REALISATIONS 

The generation of correlated Gaussian random variables is a well 
known topic. For use in the CMB this is implemented in HEALPlX 
( [Gorski et al.|2005] l, for example. Here we review the details rele- 
vant for the generation of our constrained realisations. 

In ACDM the temperature and E-type polarization are corre- 
lated as encoded in the power spectrum coefficients Cf^ , Cf^, 



and Cf^. Working in the real spherical harmonic basis we may 
generate the spherical harmonic coefficients as 



C2, 



(Al) 



(A2) 



where C,\ and C,2 are Gaussian random variables drawn from a dis- 
tribution with zero mean and unit variance. 

For our purposes we are generating constrained realisations 
of ACDM so the Ci in the above expressions are determined from 
the best fit ACDM model consistent with WMAP. Since we have 
generated co nstra ined realisations of the we have, effectively, 
chosen in | |A1| . We are left with generating the af^ . The real and 
imaginary components of it are thus generated from \A2\ using our 
known as 



C 



(CI 



<2, 



(A3) 



where ^2 is still to be chosen at random. 

APPENDIX B: DERIVATION OF STATISTICS 
FORMULAS 

BI S'^'^(0i,02) 

We wish to evaluate 5'"^'^ {61 ,62) as discussed in the text Con- 
sider the simpler case 

/X 
[c^''(e')]^d(cose) 



(Bl) 



^ (2^+l)(2f + 1) j [l-2)\{i' -2)\ ^TEr.TE 



+ 2)!(^' + 2)! 



p/(cos e)Pl (cos e) d(cos e). 



To evaluate this expression we need to perform the integral of two 
associated Legendre functions, P™ , of order m = 2, 



Pf{x)P'i,{x) Ax. 



(B2) 



Note that this integral is only defined for i, I' 2. 

For i ^ i' may proceed by directly integrating the associ- 
ated Legendre differential equation to find 

2 



l~x^ 



£{£+!) -£'(£' + 1) 



(B3) 



p!{x) 



dx 



P!,ix)qi^ 



For £ — £' more care is required. Starting from the Rodriguez 
formula 



2,d'Pi,{x) 



Pt{x)^{l-x^) 

and the recursion relation 

APi{x) ^ dP,_i(a:) 
da; Ax 
we can show that 



+ £Pi{x) 



(B4) 



(B5) 



ilf{x) = 4j)!\(a;)-4£(^-l)J,^'Va;) + ^"(^-l)^J,,4a;).(B6) 
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Here Ie,e{x) is the equivalent integral over Legendre polynomials 



encountered in the definition of Si /2 ; see Appendix A of Copi et al. 

for details. The remammg quantit.es, (.) an^lf^ 
are calculated through integration by parts and use of the recursion 
relation to find 



1 

+ 2 



+ :J{l-x[P,_i(x)]"-(2i!-l)J,_ 

and 

Jf\x) = 4-i{x) + £ [Pe^i{x)Pi{x) + 1] . 



(B7) 
(B8) 



Notice that J^^\x) is defined recursively. We can directly show 
that4^>(a;) = 0. 

With these expressions for the integrals we may write 



S^"^ {91,02) = ^c7^/,^,9(0i,e2)Cj^, 



(B9) 



This paper has been typeset from a T^U MIgX file prepared by the 
author. 



where 



If? (01,92) = 



(47r)2 
[iff, [cos 9i) 



(2^+l)(2/ + l) /(^- 2)!(^' - 2)! 



(£ + 2)!(^' + 2)! 



(BIO) 



fTQ 



(cos 02)] 



This matrix may be precomputed for rapid evaluation of 

S^'^{9i,92). 



B2 

We wish to evaluate s^'^ {9i, 02 ) as discussed in the text We 
proceed as above and consider the simpler case 

s^'^ix) = J C'^'^{9)d{cos9) (Bll) 



V (2^+1) {£-2)1 TE 



47r V {e. + 2)\ 
Pe{cos9)d{cos9). 
To evaluate this expression we need to perform the integral 
tj'^{x) = [ Pe{x)dx. (B12) 



This integral is straight forward to evaluate starting from the Ro- 
driguez formula l|B4j and integrating by parts to find 

ij'^{x) = £Pe^i{x)-{£-2)xPe{x) + 2{-l)' (B13) 

-^^[Pi+i{x)-Pe-i{x)]. 

With this we may write 

oo 

s^'^(6'i,6'2) =^Cj^iJ''(ei,6'2) (B14) 

e=2 

where 



rj''{9i,92) = "^yi^ [?r'3(coseO - l^'^{cos92)] .(B15) 

This vector may be precomputed for rapid evaluation of 

s'^'^{9i,92). 
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